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If a magnetized-orbit-charged grain encounters any abrupt inhomogeneity in plasma 
conditions during a gyro-orbit, such that the resulting in-situ equilibrium charge 
is significantly different between these regions {q\/q 2 ~ 2, where qi is the in-situ 
equilibrium charge on one side of the inhomogeneity, q 2 is the in-situ equilibrium 
charge on the other side, and q { < q 2 < 0), then the capacitive effects of charging 
and discharging of the dust grain can result in a modification to the orbit-averaged 
grain trajectory, i.e. gyro-phase drift. The special case of q\/q 2 is notioned for the 
purpose of illustrating the utility of the method. An analytical expression is derived 
for the grain velocity, assuming a capacitor approximation to the OML charging 
model. For cases in which a strong electric field suddenly appears in the wake or at 
the space-plasma-to-crater interface from solar wind and/or ultraviolet illumination 
and in which a magnetic field permeates an asteroid, comet, or moon, this model 
could contribute to the interpretation of the distribution of fields and particles. 


1. Introduction to Gyro-phase drift 

Gyro-phase drift is a guiding center drift that is directly dependent on the charging 
rate limit of dust grains (Northrop and Hill 1983). The effect of introducing a gyro- 
phase dependence on the grain charge leads to two orthogonal components of guiding 
center drift. One component, referred to here as grad-q drift (vy 9 ), results from the 
time-varying, gyro-angle-dependent in-situ equilibrium grain charge, assuming that 
the grain charging is instantaneous. For this component, the grain is assumed to 
be always in its in-situ equilibrium charge state and this state gyro-synchronously 
varies with respect to the grain’s average charge state. The other component, referred 
to here as the gyro-phase drift, arises from any non-instantaneous-charging-induced 
modification of the guiding center drift and points in the direction of — (where 
R Ld is the grain gyro-radius), i.e. the direction associated with increasing magnitude of 
in-situ equilibrium charge state. Koepke et al. (2013) simulated the effects of a gradual 
inhomogeneity on the magnitude and direction of gyro-phase drift velocity. Here we 
arrive at an analytical expression for the magnitude and direction of gyro-phase drift 
velocity for the case of abrupt inhomogeneity. 
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When there is an abrupt inhomogeneity in the plasma conditions, the in-situ 
equilibrium grain charge state will be different on either side of the inhomogeneity. If 
a magnetized-orbit dust grain gyrates through the inhomogeneity, the grain changes 
charge state q during its gyro-orbit by either gaining electrons, thus attaining a more 
negative charge state, or losing electrons, thus attaining a less negative charge state. 
The function q{t) represents the time-dependent charge of the dust grain. 

For the Orbit Motion Limited (OML) model, Goree (1994) describes a method 
for characterizing the charging time of dust grains from zero elementary charges to 
one e-folding of the in-situ equilibrium value in a homogeneous plasma as r ch = 
k T T} 12 /{ano), where T e is the electron temperature in eY, a is the dust grain radius 
in meters, no is the plasma density in m -3 , and k T , in units of seconds m -2 eV“ 2 , 
is a non-trivial function of both Ti/T e and m,/m e that is fit to a numerical model. 
It is apparent from this expression that the charging time increases with electron 
temperature, and is inversely proportional to the grain size and background plasma 
density. 

The method of Goree (1994) can be adapted to the inhomogeneous case of two 
neighboring homogeneous regions, each with its own distinct equilibrium charge state 
of the grain. Represent the two equilibrium charge states of the grain on either 
side of the abrupt inhomogeneity by qi and q 2 , with qi < q 2 < 0. The discharging 
time is the time elapsed between q(t = 0) = q\ and q(t = r d ) = q i + ^(<?i — qi), 
where ‘e’ represents the base of the natural logarithm. It is possible that q(t) will 
not reach the in-situ equilibrium charge state q 2 by the end of the discharging part 
of the gyro-orbit at time t = t d f. The charging time is the time elapsed between 
q(t d f ) = q df and q(t = t df + r c ) = q df + \{q d f — qi)- These charging time, obtained 
from OML, can fit to a discharging/charging capacitor model and the function q(t ) 
can be prescribed for the discharging/charging parts of the gyro-cycle, as shown later 
in (5.2) and (5.3). The procedure is represented graphically in Fig. 2. In principle, 
the effect of repeated excursions between abrupt inhomogeneous plasma regions is 
straightforward to analyze, although here we focus on the first excursion. 

2. Parameters, inhomogeneity, and charge model 

To demonstrate the dependence of gyro-phase drift on charging rate, an abrupt 
inhomogeneity in ultraviolet (UV) illumination is selected (UV-present for x > 0, 
UV-absent for x < 0) as the abrupt inhomogeneous condition of convenience, while 
all other plasma parameters are the same on either side of this inhomogeneity. In 
the presence of UV, a photoelectron current exists that will cause the dust grain to 
charge less negatively than in the absence of UV. Here the photoelectric current 7 P h 
is arbitrarily imposed so that q\/q 2 ~ 2. The ratio q\/q 2 ~ 2, while notional, may 
be possible using a UV source similar to Dove et al. (2012) in conjunction with the 
Magnetized Dusty Plasma Experiment (Thomas et al. 2012). 

An abrupt inhomogeneity in plasma temperature, density, or other parameters could 
also be applied in a similar fashion to this investigation, but an abrupt inhomogeneity 
in UV illumination is a readily conceivable scenario in a space or laboratory setting, 
such as the experiments of Sickafoose et al. (2000). For simplicity in the analysis, 
uniform electric and magnetic fields (E = 0 and B = Bz , with B = 4 T, respectively) 
are chosen so that only grad-q and gyro-phase drifts arise and so no other diamagnetic 
drifts occur. A magnetic field of 4 T is achievable in the Magnetized Dusty Plasma 
Experiment (Thomas et al. 2012). Isotropic UV illumination is assumed to prevent 
electron recoil and radiation forces as suggested by Rosenberg and Mendis (1995). 
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Only the motion in the x-y plane is considered, and a = 0.05 /xm grains have no 
velocity along the magnetic field direction ( v z (t ) = 0). The dust grain is started at 
the position x(t = 0) = 0, y(t = 0) = 0 so that it encounters the transition region 
(. x = 0) of this inhomogeneity at least twice during the first gyro-orbit. The dust grain 
starts with velocity components v x (t) = v_l and v y (t) = 0, where vj_ = 0.124 ms -1 
is the grain velocity perpendicular to the magnetic field. We demonstrate for the slab 

model that v± is not a crucial parameter, so v± = where k b is the Boltzmann 

constant, T N = 0.025 eV is the neutral gas temperature, and m d is the dust grain 
mass. All parameters of the problem are chosen so that V d , the dust surface potential 
relative to the local space potential or Vspace — Surface is always negative. 

For our choice of parameters in this paper, the gyro-radii of electrons and ions are 
expected to satisfy R Le , R Li > a, so the unmagnetized OML, or OML currents (Allen 
1992) are used. For Maxwellian populations of electrons and ions, and dust grain 
velocities much less than the ion thermal speed and q < 0, these currents are: 
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where e is the elementary charge, / represents the flux of photoelectrons for a given 

UV source and grain material, v the = . / AFF ; Uthi = , /AA, = n = because the 
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plasma is treated as dust-absent, V d = q/C, and C = 4ne 0 a is the grain capacitance. T e 
and f represent electron and ion temperature respectively, while m e and m, represent 
electron and ion mass respectively. For the portion of the grain’s gyro-orbit that is 
illuminated by UV (x > 0), the magnitude of the photoelectric current of electrons is 
non-zero, while it is identically zero when it is shadowed, or x < 0. The total current 
on either side of the inhomogeneity can be expressed as: 

hot = h + h + I p h for x > 0, (2.2) 


hot = h + h for x < 0. (2.3) 

Equation (2.2) describes hot when UV is present, and (2.3) describes hot when UV 
is absent. The photoelectric current in the UV-present region causes the dust grain 
to reach a less negative in-situ equilibrium charge state q 2 compared with the in-situ 
equilibrium charge state q\ in the UV-absent region. 

The photocurrent / P h is determined by 

rh 

/p h = ena 2 f = ena 2 / F (2) x W d/l, (2.4) 

Jo 

where F(l) is the flux of UV photons as a function of wavelength X for a UV source, 
x(X) is the quantum efficiency of photoemission for the grain material, and X\ is 
the wavelength cutoff, representing the longest wavelength capable of producing a 
photoelectron for a given work function (Colwell et al. 2005). In order to satisfy 
qi/qi ~ 2 for the a = 0.05-/xm grain in a plasma with background density no = 
10 13 m^ 3 and electron, ion temperatures T e = 1.6 eV, 7) = 0.025 eV, I ph ~ 2.4 x 10~ 15 
A is arbitrarily chosen. These plasma parameters are chosen because they are relevant 
to the Magnetized Dusty Plasma Experiment (Thomas et al. 2012). 
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Applying (2.2) and (2.3) to the definitions specified by (2.1) results in a single 
differential equation for q{t)\ 
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where H(x(t)) is the Heaviside step function, and x{t) is the time-dependent grain 
position along the inhomogeneity direction, x. 


3. Numerical approach to charge collection to spherical granule 

We compute the single-grain trajectory in a dust-absent plasma using a stationary 
inertial lab frame and examine the case of the abrupt inhomogeneity described in 
Sec. 2. The equations of motion for the dust grain are from the Lorentz force 
equation: 


dv x , . 

mj— = q (t) Bv y , 

(3.1) 

du y 

m d — - -q (t) Bv x . 
at 

(3.2) 


A symplectic, leapfrog integrator is used to solve the equations of motion in (3.1) 
and (3.2). There are approximately 2000 Newton timesteps (2000 A t N ) per gyro-cycle, 
with the Newton timestep held constant during the simulation. The analytical currents 
for the OML charging model (2.1) are applied to the dust grain using A q n = / tot A t c , 
where A t c is the charging timestep. This is an adaptive timestep to ensure that the 
nonlinear differential equation for the grain charge is solved correctly in accordance 
with the OML model (2.5), simultaneously with the coupled differential equations 
for the grain motion from the Lorentz force. At the beginning of the simulation, A t c 
is calculated according to A t c = e/|/ tot |, which is the time for the grain to gain or lose 
one electron. If A t c > A t N , then q(t) is held fixed at its current value and a timer 
is started with t acc = A t N . This timer, t acc , keeps track of the accumulated time that 
has passed since the last charge update. During subsequent timesteps, if t acc < A t N , 
then another Newton timestep is added to t acc . In other words, if t acc < A t c , then 
Lee = Lcc + Atjy. When the condition t acc ^ A t N is satisfied, then one electron is added 
or subtracted to q(t), depending on the sign of / tot , and the accumulated time is set 
to t acc = A t N — Remainder) At c /ALv). This way, when A t c is calculated from the new 
value of Lot at the next Newton timestep, the exact value since the last charge update 
is known. 

For this numerical scheme, q(t) is only recorded at the end of every Newton 
timestep. The total current becomes very small when the grain surface is near the 
in-situ equilibrium grain charge, and A t c becomes large. The grain fluctuates between 
neighboring charge states because, in general, C V d is not an exact integer multiple of 
the elementary charge, e. The discrete nature of the particle charging and the adaptive 
timestep are evident in the step-like plot of Fig. 3. 

The charging rate can also be arbitrarily controlled in this model by introducing 
an adjustable charging parameter, a, used in the following way: q n+ 1 = q n + A q' n , 
where A q' n = a A q as is done in Koepke et al. (2013). This is equivalent to calculating 
A t c = For a = 1, the dust grain charges without restriction, while a < 1 
implies that the grain charges more slowly. This adjustable charging rate parameter is 
discussed here only insofar as a can be used as a possible way of changing the value 
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Grain Trajectories 



Figure 1. Grain trajectories for a = 0.05 /xm. The dashed and solid lines represent the 
trajectories of P2 = T g /x ch ratios for — > oo (10 16 m -3 ) and P2 = 19.204 (for n = 10 13 m~ 3 ) 
respectively. Squares and diamonds indicate the gyro-averaged guiding centers Q = x g /x ch 
ratios for Q — > oo and £2 = 19.204 respectively. 


of x ch in simulations when other plasma parameters, such as the background plasma 
density, n 0 , are fixed. 


4. Simulation results 

Figure 1 shows a comparison of different grain trajectories for the abrupt, UV 
inhomogeneity using the OML model with instantaneous and non-instantaneous 
grain charging. Being negatively charged, the grain gyrates around B = Bz in a 
counterclockwise direction. Note how different the instantaneous-charging trajectory 
is compared with the non-instantaneous-charging trajectory. 

The dashed trajectory in Fig. 1 is the trajectory of a dust grain that experiences a 
step function inhomogeneity in UV illumination (UV light exists for x > 0, while it is 
absent for x < 0) and instantaneous charging (x g r c ), or Q — x g /x c h — > oo. For the 
portion of the grain’s gyro-orbit illuminated by UV (x > 0), there is a photoelectric 
current of electrons in addition to the plasma electron and ion currents that cause 
the grain charge to become less negative (instantaneously reaching this new charge 
state q 2 ). A smaller negative grain charge corresponds to a larger gyro-radius for 
the grain. Each returns to the shadowed region, [x < 0), cuts off the photoelectric 
current, and the grain instantaneously recovers to a more negative charge state. This 
repeating pattern of abrupt modulation of charge state results in concatenated semi- 
circles, which, in turn, give rise to a particle drift in the y-direction that is mutually 
perpendicular to the equilibrium-charge-inhomogeneity direction and the background 
magnetic field. This drift is referred to as the grad-q drift (i>y 9 ). Approximating the 
change to the gyro-radius R Ld = to occur abruptly at 0 = — rr/2, 2>xi/2, ..., which 
is a good assumption for instantaneous grain charging and x gc = x(t = 0), the grad-q 
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drift for this particle can be approximated by 


Rl , i — Rl , 2 a _ 1 v± |< 7 i | — \qi\ a 

ri/2 + T2/2' n |<7i[ + [<?2|^ 


(4.1) 


where r = 2 nm d /(qB) is the gyro-period and the subscripts 1 and 2 indicate 
equilibrium values in the UV-present and UV-absent regions respectively. Knowing 
that v± = 0.124 ms -1 , the equilibrium charge in the shadowed region is q\ = 
— 143.5 e + 0.5 e and the equilibrium charge in the UV region is q 2 = —69.5 e + 0.5 e, 
the variable charge drift is predicted to be 0.02746 ms -1 , which closely matches the 
simulation result of 0.02743 ms 1 . 

In Fig. 1, the solid trajectory is the resulting trajectory of a dust grain that has a step 
function inhomogeneity in UV illumination and non-instantaneous charging. With 
non-instantaneous charging, the dust grain does not immediately reach the new in-situ 
equilibrium charge state, and the particle trajectory is modified. In addition to the 
grad-q drift, the dust grain drifts along the inhomogeneity direction (—x -direction); 
we refer to this drift component as the gyro-phase drift. Since the gyro-center of the 
dust grain is drifting along the inhomogeneity direction, the dust grain samples the 
UV illumination region less with each gyro-orbit, causing the grain to drift out of 
the boundary associated with the inhomogeneity, and the guiding center drift velocity 
in the — x-direction decreases until it becomes zero and the gyro-center becomes 
stationary. This drift component along the inhomogeneity direction is a direct result 
of the finite charging rate of the dust grain in the model. 


5. Theoretical predictions for grad-q and gyro-phase drifts 

For the case of abrupt UV inhomogeneity and instantaneous charging, (Q — ► oo), 
the grad-q drift is given by (4.1). Since the grain instantaneously reaches the in- 
situ equilibrium grain charge, the grain velocity as a function of time for the first 
gyro-cycle can be written as 

{ v x = v± cos (co 2 t) x v y = v± sin (co 2 t) y 0 < t < ^ 

v x = v± cos x Vy = v± sin (a) it — y ^ < t < ^ ^ 

(5.1) 

assuming the grain starts with the same conditions as specified in Sec. 4. 

For the case of the abrupt UV inhomogeneity and non-instantaneous charging 
(f2 ^ oo), the particle’s velocity as a function of time in the presence of the magnetic 
field can also be computed analytically. These velocity components, being functions 
of time, can be integrated to average over a gyro-period and compute the guiding 
center drift velocity components. If the grain has an equilibrium charge of qi in the 
region where UV is absent, and the grain has an equilibrium charge of q 2 where UV is 
present, the charge on the dust grain as a function of time for the charging-discharging 
grain can be approximated by familiar equations for a discharging/charging capacitor. 
To generalize to any abrupt inhomogeneity where a grain transitions from a negative 
charge state to a less negative charge state, q\ can be associated with the charge 
of the more negative charge state and q 2 can be associated with the charge of the 
less negative charge state. For the charging part of the gyro-cycle, when the grain 
transitions to the UV-absent region, the initial charge is given by q c (t = 0) = q 2 , 
9(t = 0) = n/2 with v x (t = 0) = — v±: 

\q d (01 = (q 2 -qi)exp(-t/r d )-q 2 , 


(5.2) 
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a = 0.05-pm grain transitions from UV-absent to UV-present condition a = 0.05- pm grain transitions from UV-present to UV-absent condition 




Figure 2. Charge evolution for a = 0.05-/xm grain in the OML and capacitor models, (a) 
Charge evolution of a dust grain as it transitions from the UV-absent to UV-present region. 
The UV flux is 4 x 10 18 m _2 s _1 , which yields a photoelectron current of 5.03 x 10 _15 A. The 
capacitor model reaches \(q 2 — q i) at the same time as the OML charge model, although it 
is a different function of time, (b) Charge evolution of a dust grain as it transitions from the 
UV-present to UV-absent region. 


l<?c(0l = (9i - 92 )exp(-f/T c )-$i. (5.3) 

In the above equations, q d {t) and q c (t) are the particle charge as a function of time 
for the discharging and charging parts of the gyro-cycle. These solutions for the grain 
charge as a function of time are approximate, since the differential equation for q (t ) 
(2.5) can not be solved analytically. The actual charge evolution in the OML model 
does not precisely follow (5.2) and (5.3), but it is a reasonable approximation as shown 
in Fig. 2, which shows the grain charge as a function of time for the discharging 
and charging portions of the gyro-cycle. Also, these expressions assume that the 
grain charge is fully charged/discharged during a half gyro-cycle; this restriction can 
be lifted if q\ = q df in (5.3), where q df is the charge on the grain when the grain 
reaches x = 0, or when t = t d f. Using the following definitions: co(t) = ^ = lg ^ |g , 
the gyro-phase angle as a function of time can be calculated by integrating q d (t) and 
qft): 
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(5.5) 


(5.6) 


e c (t) = co 2 


q\ — qi 


exp — 1 + 1 


+ 0 df , 


(5.7) 


where 9 c (t = 0) = 9 df is the gyrophase angle of the grain when x = 0, after completing 
the first half of the gyro-orbit. The general 9(t) can be defined during the first 
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# <‘> = {«$*< 0- ,5 ' 8) 
With this dehnition of 9, and v(t = 0) = i>j_x, the general solution for the velocity of 
the dust grain that satisfies differential equations (3.1) and (3.2) is given by 

v x (0 = v 0x cos (9 ( t ) —6{t = 0)) - Woy sin (9 (t) - 9 (t = 0)), (5.9) 


v y (t) = vo y cos (9 ( t ) — 9 (t = 0)) + vox sin(d (f) — 9 (t = 0)), (5.10) 

where vq x and vo y are the initial velocities of the dust grain in x and y directions 
respectively. These equations can be applied to the specific conditions described in 
Sec. 2 for the charging and discharging portions of the gyro-cycle: 


v x (t) 


Vy (f) 


Vxd (*) = v± cos (9 d {t) + p) 

V xd (t Cj) COS (j? c ( t ) ®df) 

xc - -Vydit = t df ) sin (9 C (t) - 9 df ) 

v yd ( t ) = w_l cos (9 d {t) + j) 

.. = v yd (t = t df ) cos (9 C ( t ) - 9 df ) 

}L +v xd (t = t d f) sin (9 C ( t ) — 9 d f ) 


0^t ^t df 

t d f ^ t ^ t d f + t c f 


o ^ ^ t d f 

tdf ^ t ^ t d f + t c j 


(5.11) 


(5.12) 


To find the trajectory of the dust grain, the components of v must be integrated 
with respect to time. When putting (5.5) and (5.7) into (5.11) and (5.12), it is clear 
that the integrals of the form cos(h(e“* — l) + x +0) cannot be integrated analytically, 
so a numerical approach is necessary. To obtain x(t) and y(t), trapezoidal integration 
is used. Likewise, to find the guiding center drift velocity components, it is sufficient 
to numerically integrate: 


i r t,=t df i r l '=‘cf 

v xgc = ~ dt'v x ( t ') + ~ dt'v x (t'), 
tg J t’= 0 tg J ft— o 

\ r t,=t df i r t,=t cf 

VygC = — dt'Vy (V) H / dt'Uy (V), 

T g Jt '= 0 ' r g Jt '= 0 


(5.13) 

(5.14) 


where v xgc and v ygc are the x and y guiding center velocity as a function of time, t df 
is the time during the discharging part of the gyro-cycle at which the grain reaches 
the transition region (x = 0), t C f is the time at which the grain reaches 9 C = 3 / (2tc), 
and T g = t d f + t c f is the time to complete the entire gyro-orbit (the gyro-period). 
The gyro-period is defined simply as the time taken to complete one rotation in 
gyro-phase, and is not necessarily a trivial function of n/coi + n/a )2 as it was for the 
instantaneous case (5.1). Calculation of these integrals yields a guiding center drift 
of v xgc = —0.009 ms 1 and v ygc = 0.0229 ms -1 , which is close to the simulation result 
°f v xgc = —0.0088 ms -1 and v ygc = 0.0227 ms 1 . It is also clear from (5.9) and (5.10) 
that v x and v y are strictly proportional to vq x or vo y , demonstrating that the guiding 

center velocities can be normalized by v± = ^ v$ x + v% y without losing any physical 
insight. 

The simulated and analytically obtained gyro-radii ( R Ld = as a function of 

gyro-phase angle 9 are plotted in Fig. 3. As evident in these figures, the simulation 
result compares very closely with the analytical prediction based on a capacitor 
approximation of the OML charge model. 
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Figure 3. Radial distance from the origin is the gyro-radius as a function of gyro-angle. The 
grain size is a = 0.05 //m. In the analytical model, q(t ) is continuous, while ^ is discontinuous 
at 0° and ~80° because the current to the grain abruptly changes at both of these phase 
angles. 


6. Conclusions 

Gyro-motion within an abruptly inhomogeneous plasma induces charge-state 
modulation which leads to grad-q and gyro-phase drifts. Because the gyro-phase drift 
causes the grain’s gyro-center to drift along the inhomogeneity direction, dust grains 
that experience an abrupt plasma inhomogeneity with non-instantaneous charging 
eventually leave a region of structured inhomogeneity and both grad-q and gyro-phase 
drift cease. An analytical model is employed here to arrive at general and specific 
predictions of a dust grain trajectory that experiences an abrupt inhomogeneity in 
UV illumination. 

Because the gyro-phase drift causes micro-particles or nanoparticles to change 
trajectory, this effect might be used as a means of steering small grains (<100-nm 
diameter) in plasma environments where there are magnetized orbit dust grains 
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(plasma processing, tokamaks, etc.). If the spatially abrupt UV inhomogeneity is 
changed to an abrupt amplitude modulated UV inhomogeneity in time, and modulated 
at a frequency close to gyro-period of specific size grains, these grains can be moved 
in the x-y plane. Consider for example a 30-nm grain that has ~50 electrons, and a 
gyro-period of co = 48.7 rads 1 for B = 0.0875 T, which is a standard field used in 
electron cyclotron resonance (ECR) plasmas, where the electron cyclotron frequency 
is given by / ce =2.45 GHz. For a 30-nm grain population having T d = T [Ar] + =1/ 
40 eV, the thermal speed of dust grains is 0.75 ms 1 , corresponding to a gyro-radius of 
1.5 cm. If there is a UV source of sufficient intensity to cause the grains to reach the 
less negative charge state of 25 electrons, the grains will undergo a pure grad-q drift 
( r ch < r g ) if a modulation frequency of 5.8 Hz is chosen. Assuming that the phases 
of the 30-nm grains are randomized, the grains will ultimately leave the plasma 
volume with a drift direction that is perpendicular to the background magnetic field. 
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